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ABSTRACT 

Aims. The brown dwarf (BD) formation process has not yet been completely understood. To shed more light on the differences and 
similarities between star and BD formation processes, we study and compare the disk fraction among both kinds of objects over a 
large angular region in the Taurus cloud. In addition, we examine the spatial distribution of stars and BD relative to the underlying 
molecular gas 

Methods. In this paper, we present new and updated photometry data from the Infrared Array Camera (IRAC) aboard the Spitzer Space 
Telescope on 43 BDs in the Taurus cloud, and recalculate of the BD disk fraction in this region. We also useed recently available CO 
mm data to study the spatial distribution of stars and BDs relative to the cloud's molecular gas. 

Results. We find that the disk fraction among BDs in the Taurus cloud is 41 ± 12%, a value statistically consistent with the one 
among TTS (58 ± 9%). We find that BDs in transition from a state where they have a disk to a diskless state are rare, and we study 
one isolated example of a transitional disk with an inner radius of a; O.IAU (CFHT BD Tau 12, found via its relatively small mid-IR 
excess compared to most members of Taurus that have disks. We find that BDs are statistically found in regions of similar molecular 
gas surface density to those associated with stars. Furthermore, we find that the gas column density distribution is almost identical for 
stellar and substellar objects with and without disks. 
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1. Introduction 

The brown dwarf (BEQ formation process has not yet been com- 
pletely understood. Two main classes of models are currently 
envisioned to explain the origin of BDs. In the "stellar" forma- 
tion scenario (eg. Padoan & Nordlund 2002; 2004), BDs form 
like stars via gravitational collapse and fragmentation of very 
low-mass cores, followed by significant disk accretion. In the 
"ejection" formation scenario (Reipurth & Clarke 2001), BDs 
are ejected from their parent core and end up being starved for 
molecular material. In this model, accretion disks around BDs 
are radially truncated by stripping during the dynamical ejection 
(see eg. Bate 2009). There is currently no definitive evidence 
available to reject one or the other scenario either from objects' 
global properties or from study of individual objects. Studying 
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the disk proportion and structure around BDs may be a key to 
determining their formation model. 

On the one hand, Luhman et al. (2005; 2006) have stud- 
ied Spitzer Space Telescope (Werner et al. 2004) photometry of 
young objects in Taurus, and from their Spitzer 3.6 to 8 fim col- 
ors, found a BD disk fraction of x 40%, a value close to that of 
stars. Guieu et al. (2007) have studied and modeled the SEDs of 
23 substellar objects in the Taurus cloud and find that 48% of 
BDs have disks, a fraction similar to that for T Tauri stars (TTS) 
in the same region. In Chamaeleon 1, Damjanov et al. (2007) find 
no dependence of disk frequency on stellar mass in the K3-M8 
spectral type range, well into the BD regime at the age of the 
cloud (s! 2 Myr, Luhman 2004). The disk phenomenon extends 
into the planetary mass regime (M < 15 Mjup). Circumstellar 
disks have been detected around BDs within the mass range of 
giant planets (< ISMju,,), by Luhman et al. (2005, 2008a,b); 
Scholz & Jayawardhana (2008). Concerning individual objects, 
Bouy et al. (2008) have published an extensive data set on the 
substellar object 2MASS J04442713-H2512164, from the visible 
to the radio range, including the first photometric measurement 
of a BD disk at 3.7 mm, and allowing a detailed analysis of the 
disk properties. Their analysis shows that this BD has all the 
characteristics of a star with a disk. 
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On the other hand, Luhman et al. (2008a), used Spitzer pho- 
tometry on a much larger number of objects than Damjanov et 
al. (2007) in Cha I, to find a dependence on stellar mass in the 
disk fraction. In Upper Scorpius, Scholz et al. (2007) find a disk 
frequency of 37%+/-9%, a value higher than what has previously 
been derived for K0-M5 stars in the same region (at a 1 .8 cr con- 
fidence level), suggesting a mass-dependent disk frequency. 
Luhman et al. (2007) have presented a disk model for the BD 
2MASS J04381486+261 1399 showing that the disk has an inner 
hole and has an outer radius (Rom ~ 20-40 AU), a size similar to 
the predictions from hydrodynamical models (/?o„, < 10-40AU) 
by Bate et al. 2003 and Bate 2009. 

Thies and Kroupa (2007), studying recent data on the multiplic- 
ity properties of stars and BDs show them to have different bi- 
nary distribution functions. They uncovered a discontinuity of 
the multiplicity-corrected mass distribution in the very low mass 
star (VLMS) and BD mass regime, discarding a continous IMF 
with a high confidence level, and suggesting that VLMS and 
BDs on the one hand, and stars on the other, are disjoint pop- 
ulations with different dynamical histories. 
In the Taurus cloud, Guieu et al. (2007) find that BD disk frac- 
tion appears to vary with the object position within the Taurus 
filaments, unlike the known young star population. 

Disk fraction is therefore one of the keys of the very low 
mass objects formation models. With this purpose in mind, we 
present in this paper new data about BD disks in the Taurus 
cloud. We obtained new or updated Spitzer Infrared Array 
Camera (IRAC; Fazio et al. 2004) photometry data on 43 BDs 
in the Taurus cloud (Section [2] and 



3.1 1 and provide an update 



on the BD disk fraction in this region, complementing the study 
started by Guieu et al. (2007). We compare the BD disk frac- 
tion to the TTS disk fraction in the same region. We present 
a complete study of the only transition disk found among our 
BD sample (section |4]i. We also make use of recent CO data 
(Goldsmith et al. 2008) to study the distribution of stellar and 
substellar objects relative to the molecular gas (Section 5. Ill and 



the distribution of objects with and without disks (section |5T2 



2. Observations 

2.1. Brown dwarf and stellar sample 

Our BD sample results from a combination of recent studies 
dedicated to the search for BDs in the Taurus cloud. They in- 
clude studies from Briceno et al. (2002), Luhman (2004), Martin 
et al. (2002), Guieu et al. (2006), Luhman et al. (2006), and 
Luhman (2006). All BDs in our sample have been spectroscop- 
ically confirmed to be late type (> M6) and members of the 
Taurus cloud from their low surface gravity indicators and/or 
accretion lines. The stellar-substellar boundary near the spec- 
tral type M6 at the age of 1-3 Myr comes from the theoretical 
models from Baraffe et al. (1998) and Chabrier et al. (2000) and 
the current temperature / spectral type scale available for young 
low mass objects (Luhman et al. 2003). The BD surveys used 
optical and near infrared 2MASS JHR^ photometry. The com- 
pleteness of these surveys is limited by by the requirement of 
having a 2MASS counterpart, and should reach 20 - 30Mjup 
for Ay - 4 (Guieu et al. 2006; Luhman 2004). Two objects 
in the sample (J04335245-H26 12548 and J04242090-H2630511) 
have been found by their infrared excess (Luhman et al. 2006) 
using the Taurus IRAC data; we report their photometry in ta- 
ble [T] but we do not take them into account in the computation 
of the fraction of BDs with disks so that the sample is unbiased 



in terms of disks. Hence the final number of substellar members 
used in this study is 41. 

Our stellar sample includes all known Taurus members in the 
survey region at the time of this writing. 

2.2. Spitzer photometry 

Mid-infrared photometry has been obtained with the IRAC in- 
strument at 3.6, 4.5, 5.8, and 8/im on board the Spitzer Space 
Telescope. The Spitzer fluxes were extracted from the mosaic 
images obtained as part of General Observer programs 3584, 
30816 and 462 for the mapping of the Taurus cloud (Padgett et 
al. 2007; Padgett et al., in prep.). The observations were carried 
out during 2005 and 2007, with a 2x12 second high dynamic 
range frame per epoch for IRAC (two 0.4 sec exposures and two 
10.4 sec exposures). 

Data reduction and source extraction are described in detail 
in the Taurus data delivery document (Padgett et al. 2008) and 
in Padgett et al. (in prep.). A comprehensive study of the stellar 
population in Taurus using the Spitzer Space Telescope can be 
found in RebuU et al. (2009). Briefly, the IRAC basic calibrated 
data (BCDs) have been processed by the IRAC artifact mitiga- 
tion software (available on the SSC website) before being assem- 
bled in mosaics. Aperture photometry has been performed on the 
source detection output of the APEX detection algorithm. Fluxes 
have been extracted in an aperture of 2 pixel radius. The sky has 
been measured in 2 to 6 pixel radii. Fluxes were aperture cor- 
rected using the corrections listed in the IRAC Data Handbook. 
They have been converted to magnitudes given the zero magni- 
tude fluxes of 280.9 +4.1, 179.7 ±2.6, 115.0 ±1.7 and 64.1 ±0.9 
mJy for channels 1, 2, 3, and 4 respectively (Reach et al. 2005). 
The typical magnitude error is 0.05 for sources brighter than 10, 
10, 8.5, 8 mag for channels 1, 2, 3 and 4 respectively, and in- 
creases to 0.07, 0.07, 0.10 & 0.10 for the faintest sources at 14, 
13.5, 13.5, 12.7 mag. 

Guieu et al. (2007) have published Spitzer photometry for 
23 BD in the Taurus cloud based on the original Spitzer Taurus 
survey. Since then, the Spitzer survey has been extended to a 
larger region (Padgett et al. 2009) and new data are now avail- 
able for more BDs. In Table[T| we list 43 BDs currently known in 
the Taurus cloud with Spitzer IRAC photometry. The data from 
Guieu et al. (2007) have been (slightly) updated and are there- 
fore superseded by these new values. Among the 20 additional 
objects reported here, listed in table[r| 9 show an infrared excess. 



taken as a signature of a circumstellar disk (see section 3.1 1. The 



table also includes visual absorption values (col 10). These val- 
ues come from the literature when available, they are often con- 
verted from Ay to Av using Aj - 0.265 Ay, or estimated from 
their spectral types combined with their J-H color 



2.3. CO data 

In order to study the distribution of BD with and without disks 
with respect to the surrounding molecular cloud, we use the re- 
cently published CO data from Goldsmith et al. (2008). They 
have used '^CO and '^CO data to compute the molecular hy- 
drogen column density, N(H2), over a large area of the Taurus 
cloud, encompassing all the BD of our sample. The distribution 
of N(H2) at the position of every BD of the sample is listed in 
table [T] (see section 



5.1 for a description of how the density is 



computed for every object). 
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Table 1. IRAC Spitzer photometry for BDs. " denotes the two objects found from their IR excess. 



Name 


2MASS name 


N(H2) 


Sp Type 


Disk 


[3.6] 


[4.5] 


[5.8] 


[8] 


Av 






(10^' cm-2) 




presence 












J041411+281153 


04141188+2811535 


5.89 


M6.25 


yes 


10.93 


10.40 


10.12 


8.99 


1.06 


KPNO-Tau 1 


04151471+2800096 


3.60 


M9.00 


no 


13.23 


12.72 


12.94 


12.81 


0.41 


KPNO-Tau 2 


04185115+2814332 


5.88 


M6.75 


no 


12.25 


12.11 


12.02 


11.84 


0.37 


KPNO-Tau 12 


04190126+2802487 


4.98 


M9.00 


yes 


13.97 


13.61 


13.23 


12.75 


0.53 


CFHT-Tau 14 


04221644+2549118 


1.86 


M7.75 


no 


11.48 


11.34 


11.28 


11.23 


0.57 


CFHT-Tau 9 


04242646+2649503 


2.04 


M6.25 


yes 


11.16 


10.88 


10.51 


9.83 


0.91 


J042728-^261205 


04272799+2612052 


2.03 


M9.50 


no 


12.57 


12.37 


12.21 


12.08 


2.45 


CFHT-Tau 15 


04274538+2357243 


1.88 


M8.25 


no 


13.24 


13.15 


13.25 


13.06 


1.30 


KPNO-Tau 5 


04294568+2630468 


1.92 


M7.50 


no 


11.05 


11.02 


10.94 


10.83 


0.00 


KPNO-Tau 6 


04300724+2608207 


2.39 


M9.00 


yes 


13.12 


12.77 


12.42 


11.58 


0.88 


CFHT-Tau 16 


04302365+2359129 


2.42 


M8.50 


no 


13.23 


13.16 


13.04 


12.99 


1.51 


KPNO-Tau 7 


04305718+2556394 


2.25 


M8.25 


yes 


12.62 


12.29 


11.99 


11.25 


0.00 


CFHT-Tau 13 


04312669+2703188 


3.83 


M7.25 


no 


12.90 


12.75 


12.72 


12.67 


3.49 


CFHT-Tau 7 


04321786+2422149 


5.43 


M6.50 


no 


9.98 


9.87 


9.76 


9.72 


0.00 


CFHT-Tau 5 


04325026+2422115 


6.10 


M7.50 


no 


10.46 


10.27 


10.10 


10.07 


9.22 


CFHT-Tau 12 


04330945+2246487 


3.72 


M6.50 


yes 


10.86 


10.63 


10.34 


9.95 


3.44 


CFHT-Tau 1 


04341527+2250309 


3.56 


M7.00 


no 


11.24 


11.10 


10.98 


11.02 


3.10 


CFHT-Tau 1 1 


04350850+2311398 


2.16 


M6.75 


no 


11.19 


11.12 


11.04 


10.99 


0.00 


KPNO-Tau 9 


04355143+2249119 


2.03 


M8.50 


no 


13.63 


13.52 


13.67 


13.41 


0.00 


CFHT-Tau 2 


04361038+2259560 


3.07 


M7.50 


no 


11.63 


11.43 


11.34 


11.33 


0.00 


CFHT-Tau 3 


04363893+2258119 


2.55 


M7.75 


no 


11.79 


11.69 


11.60 


11.57 


0.00 


1043801-^255857 


04380083+2558572 


4.40 


M7.25 


no 


9.60 


9.47 


9.37 


9.31 


0.64 


J043814+261139 


04381486+2611399 


4.20 


M7.25 


yes 


10.80 


10.21 


9.64 


8.92 


0.00 


GMTau 


04382134+2609137 


5.24 


M6.50 


yes 


9.27 


8.77 


8.43 


7.81 


4.34 


J043904+254426 


04390396+2544264 


4.88 


M7.25 


yes 


10.75 


10.46 


10.02 


9.14 


0.40 


CFHT-Tau 4 


04394748+2601407 


6.21 


M7.00 


yes 


9.54 


9.07 


8.60 


7.79 


2.64 


J0441 11+255511 


04411078+2555116 


5.07 


M6.50 


yes 


10.78 


10.35 


9.92 


9.22 


1.80 


J044148+253430 


04414825+2534304 


5.84 


M7.75 


yes 


11.43 


10.93 


10.50 


9.54 


1.06 


J044427+251216 


04442713+2512164 


3.10 


M7.25 


yes 


9.56 


9.00 


8.36 


7.43 


0.00 


J041524+291043 


04152409+2910434 


2.26 


M7.00 


no 


11.86 


11.79 


11.66 


11.49 


1.18 


J041618+275215 


04161885+2752155 


4.83 


M6.25 


no 


10.88 


10.78 


10.67 


10.68 


1.07 


J042900+275503 


04290068+2755033 


1.99 


M8.25 


yes 


12.30 


11.99 


11.60 


10.92 


0.00 


J043 119+233504 


04311907+2335047 


2.27 


M7.75 


no 


11.66 


11.53 


11.56 


11.46 


0.90 


J043203+252807 


04320329+2528078 


1.91 


M6.25 


no 


10.30 


10.20 


10.13 


10.09 


0.00 


J043223+240301 


04322329+2403013 


3.01 


M7.75 


no 


10.89 


10.83 


10.79 


10.68 


0.00 


JU4j j4z+z5/o4 / 


043 J4z9 1 +25264 /O 


2.71 


M8.75 


no 


12.76 


12.63 


12.52 


12.47 


0.62 


J043545+273713 


04354526+2737130 


1.65 


M9.25 


no 


13.18 


13.11 


12.93 


13.07 


0.28 


J043610+215936 


04361030+2159364 


1.08 


M8.50 


yes 


13.02 


12.74 


12.41 


11.74 


0.33 


1042420+263051° 


04242090+2630511 


3.77 


M6.50 


yes 


11.83 


11.45 


11.00 


10.40 


0.85 


1043352+261254" 


04335245+2612548 


3.78 


M8.50 


yes 


13.17 


12.66 


12.28 


11.40 


5.20 


J042630+244355 


04263055+2443558 


3.37 


M8.75 


yes 


12.57 


12.21 


11.76 


11.08 


0.00 


J042 154+265231 


04215450+2652315 


5.30 


M8.50 


no 


13.22 


13.12 


12.90 


12.80 


2.97 


CFHT-Tau 10 


04214631+2659296 


3.53 


M6.25 


yes 


11.54 


11.32 


11.05 


10.45 


3.59 



3. Brown dwarf disks 

3. 1 . Update of the brown dwarf disk fraction 

We use the Spitzer color information to decide whether a BD (or 
a TTS) has a disk or not. Following Gutermuth et al. (2008), all 
Taurus members with the following color constraints are likely 
BDs or stars with disk. 

C [4.5] - [8] > 0.5 

& 

[3.6] - [5.8] > 0.35 (1) 
& 

, [3.6] - [5.8] < ^ X ([4.5] - [8] - 0.5) + 0.5 

In fig[T] we have plotted the [4.5]-[8] vs. [3.6]-[5.8] Spitzer 
color-color data for all the BD from table [T] (triangles) superim- 
posed on the same plot for previously known Taurus stars with 
(x symbols) and without (+ symbols) disks. The grey area delin- 
eates the objects with a disk. Apart from some objects at the edge 



of this zone, objects (TTS and BDs) with and without infrared 
excess are well separated in the plot. 

Almost all the BDs in table[T]have Ay < 4 (one has Ay - 5.2 
and one has Ay = 9.2). We choose to restrict our analysis to the 
part of the cloud with Ay < 4, for BDs and for stars. Removing 
the two objects found from their disk excess and mentioned in 
section [2T| and two more from this absorption constraint yields 
a BD disk fraction of 16/39 = 41 ± 12 %. If we use the Spitzer 
photometry to compute the fraction of disks in the TTS, we find 
that among 103 stars with Ay < 4, 60 have disks, yielding a disk 
fraction of 58 + 9 %. The disk fraction of TTS and BD in Taurus 
thus appears identical at the one sigma level. 

3.2. Brown dwarf transition disks 

The plot in figure [T] shows two well separated populations. 
However, the plot is crowded and makes the reading of possi- 
ble transition objects difficult. In order to ease the measurement 
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[3.6]-[5.8] 

Fig. 1. Taurus TTS and BDs in a Spitzer IRAC color-color dia- 
gram. Stars without disks are plotted as 'x' symbols while stars 
with disks are plotted as '+' symbols. BDs are plotted as empty 
or solid triangles depending on their colors / disk status. 



of the color excess, we follow Damjanov et al (2007) in plotting 
[4.5]-[8] vs. J-H, this latter color being used as a proxy for the 
objects photosphere, for all the sources of our sample in figure|2] 
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Fig. 2. Taurus color-color diagram in the [4.5] - [8] /7 - // plane. 
Symbols are identical to those in figure [T] Note that this figure 
is focused on the J-H BD regime. Twenty more TTS are dis- 
tributed between J - H-2 and / - H-3.5 



While the values of J-H span a wide range of values due to 
the distribution of spectral types, the [4.5] -[8] colors fall into 
two groups, depending on whether or not the objects have disks. 
The objects with disks show more scatter in their [4.5] -[8] index 
since it is dominated by disk excess emission. Figure |2] shows 
some transition objects and confirms that CFHT-Tau 12 is a BD 
with a disk of this type. We study this peculiar object in more 
details in the following section. 



4. CFHTTau 12 

CFHT Tau 12 presents a very peculiar SED, with an excess 
starting at 8/im and extending to longer wavelengths. In or- 
der to ascertain the nature of the circumstellar environment of 
CFHT Tau 12, we have performed a detailed modeling of its 
SED. 



4.1. Disk model 

We use the 3D Monte-Carlo continuum radiative transfer code 
MCFOST (Pinte et al. 2006). Our model includes multiple scat- 
tering, dust heating assuming radiative equilibrium, and contin- 
uum thermal re-emission. 

We use a density distribution with a Gaussian vertical profile 
p{r,z) - Pair) exp(-z^/2 /i^(r)), assuming a vertically isother- 
mal, hydrostatic, non self-gravitating disk. Power-law distribu- 
tions are utilized for the surface density E(r) - I.{)(r/ro)"' and 
the scale height h(r) = hoif/ro'f where r is the radial coor- 
dinate in the equatorial plane, ho the scale height at the radius 
ro - 100 AU. The disk extends from an inner radius to 
an outer limit radius rout = 300 AU. The central star is repre- 
sented by a uniformly radiating sphere with a NextGen spec- 
trum from Allard et al. (2000) (reft = 2 900K, logg = 3.5) and 
an Ay = 3.44. 

We consider homogeneous spherical grains, using the dielec- 
tric constants described by Mathis & Whiffen (1989) in their 
model A. The differential grain size distribution is given by 
dn(fl) oc a^^'^ da with grain sizes between flmin = 0.03//m and 
flmax, which is taken as a free parameter. The mean grain density 
is 0.5 g.cm"^ to account for fluffiness. Extinction and scatter- 
ing opacities, scattering phase functions and Mueller matrices 
are calculated using Mie theory. Dust and gas are assumed to be 
perfectly mixed, and grain properties are taken to be independent 
of position within the disk. 

4.2. Model fitting 

Due to the ambiguities with SED fitting, it is not possible to 
constrain model parameters independently. The robust estima- 
tion for the range of validity of the parameters instead requires 
the potential correlations between each of the parameters to be 
taken into consideration. With this in mind, we systematically 
explored a grid of models by varying 7 free parameters whose 
values are listed in Table |2] 

Comparisons between the models and the observations were 
drawn according to reduced calculations. To determine the 
range of validity for the parameters, we used a Bayesian infer- 
ence method (Press et al. 1992; Lay et al. 1997; Pinte et al. 2007; 
Pinte et al. 2008). This technique allows us to estimate the prob- 
ability of occurrence of each parameter value. The relative prob- 
ability of a single point of the parameter space (i.e. one model) 
is proportional to exp{-x^ /2), where refers to the reduced 
of the corresponding model. All probabilities are normalized at 
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Table 2. Range of values explored for each parameter 



Parameter 


min value 


max value 


# values 


sampling 


n„ (AU) 


0.01 


100 


20 


log 


Md„s, (Mo) 


10-'^^ 


lO"*^ 


8 


log 


flmax (AU) 


1 


1000 


10 


log 


P 


1.0 


1.25 


5 


linear 


a 


-0.5 


-1.5 


5 


linear 




10 


30 


10 


linear 


cos(i) 


0.05 


0.95 


10 


linear 



the end of the procedure so that the sum of the probabilities of 
all models over the entire grid is equal to 1 . 

The Bayesian method relies on a priori probabilities for the 
parameters. Here, we assume that we do not have any prelim- 
inary available information, and we choose uniform a priori 
probabilities which correspond to a uniform sampling of the pa- 
rameters. However, in the absence of any data, some grid points 
are more likely than others: consideration on solid angles show 
that an inclination between 80 and 90° (close to edge-on) is 
more likely than a inclination between and 10°. Uniformly 
distributed disk inclinations and orientations in three dimensions 
correspond to a uniform distribution in the cosine of the inclina- 
tion. Also, some physical quantities, like the inner radius, tend 
to be distributed logarithmically. The grid was built according to 
these distributions (Table|2|. 

4.3. Results 

Figure [3] presents the best fit SED (i.e. the model with the 
smallest ;^'^) and the relative figures of merit estimated from the 
Bayesian inference method for the disk inner radius, and for the 
inner radius and maximum grain size simultaneously. These re- 
sults were obtained from marginalization (i.e. summing) of the 
probabilities of all models, where one (or two) parameter is fixed 
successively to its different values. The resulting histograms in- 
dicate the probability that a parameter takes a certain value, 
given the data and assumptions of our modeling. The width of 
the probability curve in the central panel is a strong indicator of 
how well the disk mass is constrained. 

The range of validity for the disk mass is defined as the in- 
terval [Ml, M2] where: 

p(r)dr^ p(r)dr=—^ (2) 

JM„i„ Jm2 ^ 

with y = 0.997. The interval is [5 lO^'^ 5 10"'' Mq] (Figjs] cen- 
tral panel) and corresponds to a 99.7 % confidence interval, i.e. 
equivalent to a 3 cr interval in the case of a Gaussian distribu- 
tion of probability. This constraint is set by the upper limits at 
70 and 160yL(m. Unsurprisingly, most of the parameters are not, 
or only slightly, constrained. For instance, the same analysis for 
the inner radius (Fig [3] right panel) gives an interval of [0.01, 
1.5] AU. 

CFHT Tau 12 is surrounded by a very low mass disk, with 
a dust mass lower than 1.510^3 Earth mass. This is signifi- 
cantly lower than disk masses measured so far for BDs (Scholtz 
et al. 2006). The mass we derive is comparable to masses of 
debris disks around M stars, which have been estimated to be 
between 10"^ and few earth masses (see for instance Fig. 3 of 
Wyatt 2008). Because the disk excess is only detected in the 
mid-infrared, an hypothesis is that most of the outer disk has 
been removed (potentially via an ejection mechanism of the BD 
for instance). 



5. Spatial distribution of brown dwarfs in the Taurus 
cloud 

5.7. Stars and brown dwarfs relative to Taurus molecular gas 

Luhman (2006) studied the nearest neighbor distance and found 
no differences between the star and BD distributions. Such a re- 
sult is consistent with recent computations by Bate (2009) who 
finds no significant mass segregation. However, models of spa- 
tial distribution of BDs produced by the decay of small-N stellar 
systems by Goodwin et al. (2005) show that similar distributions 
sometimes arise even with ejections. In other words, a different 
spatial distribution for stars and BDs is probably a signature of 
ejections, but a lack of difference does not necessarily exclude 
the ejection scenario. 

As the star density distribution is only a proxy for the actual 
density structure of the cloud where the objects were born, and 
in order to better track the relation between stars & BDs and the 
underlying molecular gas, we have used the high spatial resolu- 
tion molecular H2 surface density map computed by Goldsmith 
et al. (2008) from recent '^CO and '^CO data, and we study the 
distribution of stellar and substellar objects relative to the actual 
cloud. We stress that strictly speaking, the surface density (cr) is 
only a proxy for the volume density (p). However, a given ob- 
ject cannot be close to a volume density peak while far from a 
surface density peak. We measure the value of the H2 column 
density at the position of every object, star or BD, in the Taurus 
cloud, limiting ourselves to objects with Ay < 4 as discussed 
above. Each measurement encompasses x 3000 x 3000 AU at 
the Taurus distance. Figure |4] shows the histogram of the un- 
derlying H2 column density for stars (solid line) and BDs (dot- 
dashed line). The BD distribution shows a sharper peak at low 
surface densities, but overall the two histograms are very simi- 
lar. We have performed a KS test on the two distributions and 
find that there is 36% probability that the two samples are drawn 
from the same distribution. This result is consistent with previ- 
ous studies of the nearest neighbor distance and shows that the 
distributions of both kinds of objects are the same with respect 
to the underlying molecular cloud. 

5.2. Brown dwarfs with and without disks with respect to 
Taurus molecular gas 

In this section, we use N(H2), the H2 surface density data to 
study the distribution of objects with and without disks relative 
to the molecular cloud. We apply the method to stars and BDs. 
Such a study is motivated by the fact that if TTS with and with- 
out disks appear relatively regularly distributed across the Taurus 
filament. However, there is a possibility that BDs with disks are 
more frequent in the northern part of the cloud (Guieu et al., 
2007). If we split the Taurus filaments into two parts containing 
approximately the same number of objects (see dashed line in 
figure |6]l, and compute the disk fractions for TTS and BD in the 
northern and the southern filaments, we find that, while the TTS 
disk fraction is almost identical in both regions, the northern BD 
disk fraction is significantly higher than the southern one (see 
Table[3]l. 

As in section |5.1| we compute the underlying H2 column 
density at the position of every object in the cloud (stars and 
BDs), with and without disks. The corresponding distribution of 
N(H2) for BDs (with and without disk) and stars (with and with- 
out disk) is shown in figure |5] We then perform a KS test to 
estimate the probability that for the two kinds of objects (stars + 
BDs), both populations (with and without disks) are drawn from 
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Fig. 3. SED modeling of CFHT Tau 12. Left: SED and best model. Crosses: data points, triangles: upper limits. Full blue line: 
synthetic SED. Red line: stellar photosphere. Center: marginal Bayesian probability for the disk mass. Right: marginal Bayesian 
probabihty for the inner radius. 



0.4 

i 

;-i 

2 0.3 



s 

o 



0.2 



0.1 



0.0 



r— |- 



: \ I 



_.+._: 



i 
i 



N(H2) (1021 



Fig. 4. N(H2) histogram for BDs and TTS with Ay < 4. SoUd 
line : stars, dot-dashed line : BDs; the vertical lines indicate the 
+ 1 cr error bars. 



the same distribution. We find that the probability that each kind 
of object (BDs or stars), with and without disks are drawn from 
the same distribution is very high (18% for BDs, 31% for stars). 
The distribution of objects with or without disks appears inde- 
pendent of the underlying gas distribution, suggesting that the 
suiTounding molecular cloud physical properties (at least den- 
sity) are not directly linked to the existence of a disk around a 
newly formed object. Thus if the peculiar BD distribution men- 
tioned by Guieu et al (2007) is real, it is not linked to the un- 
derlying cloud gas density. In a study on WTTS and CTTS in 
Taurus, Bertout et al. (2007) conclude that the disk lifetime is 
4x10* (M«/Mq)"^^ yr in this cloud. A variation of the disk life- 
time with the central object' mass associated with small age 
difference between the filaments would be a possible explana- 
tion for the disk fraction difference between the northern and the 
southern part of the Taurus cloud. However, a shorter disk life- 
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Fig. 5. Molecular hydrogen surface density distribution for stars 
(upper panel) and BDs (lower panel); in each plot, the solid line 
is the distribution of object with disk, and the dashed line, with- 
out disk. 



time for BDs is very difficult to reconcile with BD disk fraction 
measurements in older star forming regions like Cha I, IC 348 
and cr Ori (Luhman et al. 2005; 2008a,b) where it is similar to 
that of stars. One possible explanation for the different disk frac- 
tions among stars and BDs between the north and south (if it is 
real) is that the southern clouds are indeed older, but the disk 
fraction for stars is artificially high due to another explanation. 
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Fig. 6. Taurus TTS and BDs superimposed on CO molecular gas. 
The dashed line delineates the separation chosen between the 
northern and southern filaments (see text). Symbols are identical 
to those in figure [T] 



Surveys of Taurus have been biased toward finding stars with 
disks rather than without disks. Many of the previous wide-field 
surveys have selected members based on disk indicators like 
Ha and IR excess emission, whereas objects without disks are 
found by proper motions and the all-sky ROSAT survey, which 
are rather shallow. In comparison, the BDs have been identified 
mostly through optical color-magnitude diagrams, which are un- 
biased in terms of disk presence. 



Table 3. CTTS and BD disk fraction in the northern and southern 
Taurus filament for objects with Ay < 4. 



Filament 


TTS 




BD 




North 


37/68 


= 54 ± 9 % 


14/25 = 56 


± 19% 


South 


23/35 


= 66 ± 14 % 


3/14 = 21 ± 


14% 



Testing the distribution of objects in the cloud relative to the 
gas was motivated by the fact that previous papers like Goodwin 
et al. (2005) obtained similar nearest neighbor distributions, even 
with ejections. The main result from our comparison of the BD 
and star distributions with respect to the underlying molecular 
gas is that BDs and stars are most probably drawn from the 
same distribution (p value of the null hypothesis is 0.36). This 
work complements and yields a picture consistent with the one 
published previously by Luhman (2006) in which BD and stars 
are distributed in a similar manner across the Taurus cloud, and 
present similar nearest neighbor distributions. The presence of » 
two X more BD than stars in the first 2 x 10^' cm"^ H2 density 
bin remains to be explained. Yet, as both distributions appears 
to be statistically identical, this might be due to small number 
statistics effect. 



6. Conclusions 

We have presented new photometric measurements in the 3 - 
8yum domain from Spitzer IRAC instrument, complementing 
previous IR photometry by Guieu et al. (2007) for more than 40 
BDs in the Taurus cloud. We associate these IR measurements 
with new H2 surface density computations from Goldsmith et al. 
(2008) to study and compare the distribution of stellar and sub- 
stellar objects relative to the underlying molecular gas. Based on 
this larger sample, we update the disk fraction for stars and BDs 
in the Taurus cloud and find that the disk fraction around BDs is 
similar to the disk fraction in their TTS counterparts. 

We find that transition disks are very rare among BDs, a 
common result with TTS showing that the disk's disappearance 
timescale is also very short for substellar objects. These results 
call for further theoretical models of dust coagulation and set- 
tling as a function of the central object's mass. We model the 
only transition disk we find around CFHT BD Tau 12, an ob- 
ject that could have experienced more dust evolution than other 
members in Taurus. Note that Luhman et al (2006) find a slightly 
earlier spectral type than us for this object (M6 instead of M6.5) 
so CFHT BD Tau 12 could be an object at the border between 
VLMS and BD rather than a bona fide young BD. 

We study the distribution of stars and BDs relative to the H2 
molecular gas surface density computed from CO observations 
and find that both distributions are likely drawn from the same 
distribution, even if there appears to be an excess of BDs found 
at low H2 surface density. A similar study for objects with and 
without disks shows that they have the same distribution relative 
to the underlying gas. 

Finally, the finding that the fact that BD with disks are less 
numerous in the southern than in the northern portion of the 
Taurus filaments could be a reflection of an older age for the 
southern part of the cloud combined with a shorter lifetime for 
BD disks. A further determination of the ages of the Taurus fil- 
aments via HR diagram for instance is difficult due to the large 
dispersion one obtains in the ages. Disk fraction measurements 
in Cha I or IC 348 show no evidence that the lifetimes of sub- 
stellar objects are shorter 
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